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ABSTRACT 

The strategy followed so far in the performed or proposed tests of the gen- 
eral relativistic Lense-Thirring effect in the gravitational field of the Earth with 
laser-ranged satellites of LAGEOS type relies upon the cancelation of the dis- 
turbing huge precessions induced by the first even zonal harmonic coefficient J2 
of the multipolar expansion of the Newtonian part of the terrestrial gravitational 
potential by means of suitably designed linear combinations of the nodes Q of 
more than one spacecraft. Actually, such a removal does depend on the accu- 
racy with which the coefficients of the combinations adopted can be realistically 
known. Uncertainties of the order of 2 cm in the semimajor axes a and 0.5 
milliarcseconds in the inclinations / of LAGEOS and LAGEOS II, entering the 
expression of the coefficient Ci of the combination of their nodes used so far, 
yield an uncertainty 6ci = 1.30 x 10~^. It gives an imperfectly canceled J2 sig- 
nal of 10.8 milliarcseconds per year corresponding to 23% of the Lense-Thirring 
signature. Uncertainties of the order of 10 — 30 microarcseconds in the inclina- 
tions yield 6ci = 7.9 x 10~^ which corresponds to an uncanceled J2 signature 
of 6.5 milliarcseconds per year, i.e. 14% of the Lense-Thirring signal. Concern- 
ing a future LAGEOS-LAGEOS II-LARES combination with coefficients ki and 
^2, the same uncertainties in a and the less accurate uncertainties in I as be- 
fore yield 6ki = 1.1 x 10~^,6k2 = 2 x 10~^; they imply a residual J2 combined 
precession of 14.7 milliarcseconds per year corresponding to 29% of the Lense- 
Thirring trend. Uncertainties in the inclinations at ^ 10 microarcseconds level 
give Ski = 5 X 10~^, 6k2 = 2 x 10~^; the uncanceled J2 effect is 7.9 milharcseconds 
per year, i.e. 16% of the relativistic effect. 

Subject headings: Experimental tests of gravitational theories — Satellite orbits — 
Harmonics of the gravity potential field; geopotential theory and determination 
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PACS: 04.80. Cc, 91.10.Sp, 91.10.Qm 



1. Introduction 



According to tlie Einsteinian General Theory of Relativity (GTR), the Lense-Thirringll 
precession of the longitude of the ascending nodcl fi of a test particle orbiting a cent ral, 
slowly rotating body of mass M and angular momentum S is fJLense fc ThirringlllQlSl ) 



VL 



2GS 



LT 



c2a3(i_e2)3/2' 



where G is the Newtonian constant of gravitation, c is the speed of light in vacuum, a, e are 
the semimajor axis and the eccentricity, respectively, of the test particle's orbit; note that 
eq. ([I]) does not depend on the inclination / of the orbit to the central body's equator. The 
Lense-Thirring effect is a consequence of the fact that, in its weak-field and slow- motion 
approximation, the Einstein's field equations of GTR get linearized, thus resembling the 
linear equations of the Maxwellian electromagnetism. In such a framework, analogously to 
the magnetic field induced by electric currents acting on a moving electric charge through 
the Lorentz force, mass-energy currents like those of an isolated rotating mass generate a 
gravitomagnetic field in the space surrounding it which acts on a moving test particle with 
a non-central, Lorentz-like force perturbing its Keplerian motion. 

Attempts to detect the Lense-Thirring node precessions in the gravitational field of 
the Earth have been performed so far with the passive geodetic LAGEOS and LAGEOS II 



satellites (ICiufo 



(Pearlman et al. 



ini et al.ll2009l ) tracked with the Satellite Laser Ranging (SLR) technique! 



20021) : a total accuracy of the order of approximately 10% is claimed 
f lCiufolini et al.l|2009h . A further LAGEOS-li ke SLR target name d LARES, should be 
launched in latqj 2011 with a VEGA rocket ( ICiufolini et al.ll2009l ): its proponent claims 



that, in conjunction with the existing LAGEOS and LAGEOS II, it will be able to reach a 



^According to iPfisterl (120071 ) . it would be more correct to speak about an Einstein- 
Thirring-Lense effect. 

^It is an angle in the {xy} reference plane, coinciding with the equatorial plane of the 
central body, between the reference x axis and the line of the nodes, which is the intersection 
between the test particle's orbital plane and the reference plane. 

■^It allows to precisely measure the range p between a laser station and a spacecraft that 
is equipped with retrorefiectors hke just the LAGEOS satellites. The range is deduced from 
the elapsed time of flight for a pulse of laser light traveling from the ground station to the 
satellite and back again. 

■^See on the WEB: http://spaceflightnow.com/tracking/index.html 
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Table 1: Orbital parameters and Lense-Thir ring; node precessions of L AGEOS, LAGEOS II 



and LARES for 5® = 5.86 x lO^^ kg m^ s"^ flMc Carthv fc Petitll2004f ). The semimajor axis 



a is in km, the inclination I is in deg, and the Lense-Thirring rate ^lt is in nias yr 



Satellite 



Q 



LT 



LAGEOS 
LAGEOS II 
LARES 



12270 
12163 

7828 



0.0045 

0.014 

0.0 



109.9 
52.65 
71.5 



30.7 
31.5 
118.1 



~ 1% accuracy in measuring the Lense-Thirring effect fjCiufolini et al.ll2009l ). The values of 
the relevant orbital parameters and of the Lense-Thirring node precessions for LAGEOS, 
LAGEOS II and LARES are in Tabled! They amount to a few ten-hundred milliarcseconds 
per year (mas yr~^ in the following) corresponding to linear shifts of about 2 — 5 m per year 
at their altitudes. 

Actually, the nodes of such satellites are affected by much larger secular precessions 
CLj^ = CliJi induced by the even {i = 2,4,6,...) zonal (m = 0) harmonic coefficients 
Je,£ = 2,4,6,... of the multipolar expansion of the Newtonian part of the terrestrial 
gravitational potential which account for the departures from spherical symmetry of the 
Earth because of its diurnal rotation (JTapley et al.ll2004f ). The even zonal harmonics, 
defined as Ji = — \/2i + 1 Cm,i = 2, 4, 6, ... in terms of the normalized Stokes coefficients 
Ceo,(^ = 2,4,6,... (JTapley et al.l 12004 ). are directly estimated as solve-for parameters in 
global Earth's gravity field solutions f obtained by processing huge data sets from dedicated 
satellite-based mission like CHAMFO and, especially, GRACEjU. The most effective even 
zonals in perturbing the satellite s' nodes are the low-degree ones; the coefficients Q,e of the 
node precessions for i = 2,4 are (jCiufolinilll996l : Iloridl2003l ) 






Qa = Q- 



(f) 



(l-e2)2 



2 1+1 



(l-e2) 



^lyl 



(7sin2/-4) 



(2) 



where n = ^jGMja^ is the satellite's Keplerian mean motion and R is the mean equatorial 
radius of the central body; contrary to the Lense-Thirring precession of eq. (fTj), the classical 
precessions of eq. ([2]) depend on the inclination J. The node precessions due to the first 
even zonal J^ for LAGEOS, LAGEOS II and LARES are listed in Table El It can be noted 
that they are 7 orders of magnitude larger than the Lense-Thirring precessions of Table [H 



^They are publicly available on the WEB at http://icgem.gfz-potsdam.de/ICGEM/ 

^See on the WEB: http://www-app2.gfz-potsdam.de/pbl/op/champ/ 

'^See on the WEB: http://www-app2.gfz-potsdam.de/pbl/op/grace/index_GRACE.html 
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Table 2: Node precessions Qj^ = ^.2^/2, in nias yr ^, of LAGEOS, LAGEOS II and 
LARES due to J2 We used C20 = -4.841692151273 x 10 "^ from the ITG-Grace2010s 
f lMaver-Giirr et al.lboioh global solution. Recall that Je = -V^ITT Cio, £ = 2, 4, 6, .... 



LAGEOS 



LAGEOS II 



LARES 



n 



J2 (mas yr" 



4.516313623 x 10" 



-8.303250890 x 10* 



-2.0298203351 x 10^ 



Thus, suitable linear combinations of the nodes of more than one satellite have been set 
up in order to purposely cancel out, by construction , the impact of one or more even zonals 
according to a strategy put forth by ICiufolinil (119961 ). In particular, the tests pe rforrned so 



far have been conducted with the following LAGEOS-LAGEOS II combinatioiij (jl 

/•(2L) ^qLAGEOS _|_ qLAGEOS II 



Lorio 



20061 



where 



Cl 



^LAGEOS 
'^LAGEOS II 



COS/i 



LAGEOS 



COS /lAGEOS II 



^LAGEOS 



^LAGEOS 



")■( 



QlAGEOS II 
^LAGEOS 



7/2 



(3) 



(4) 



The f uture combi nation involving LARES as well, designed to remove the effect of J2 and 
J4, is (jlorioll2005[ ) 

J(3L) ^^LAGEOS ^ ^^^LAGEOS II ^ ^^f^^^^^^S^ (5) 

in which 



oLARESnLAGEOS oLAGEOSfSLARES 

t, ".2 ".4 ^'.2 ".4 

"'I qLAGEOS IIj^LARES_QLARESf2LAGEOS II) 



fSLAGEOSoLAGEOS II oLAGEOS IIALAGEOS 
"'2 j^^LAGEOS IIj^LARES_QLARESf2LAGEOS II • 



(6) 



It is analogous to the cor nbination of th e nodes of LAGEOS and LAGEOS II and the 
perigee u of LAGEOS II f lCiufolinilll996l ) used in the earlier tests fICiufolini et al.lll998h : 
the coefficients of the precessions of the perigee of LAGEOS II have to be replaced by those 
of the precessions of the node of LARES. Tabled! eq. (|2]), eq. @, and eq. (Q yield the 
numerical values of Ci,ki, k-2, shown in Table El Table |31 eq. (E]) and eq. ^ yield 47.8 mas 
yr~^ and 50.8 mas yr~^, respectively, for the predicted Lense-Thirring combined precessions. 
The LAGEOS-LAGEOS II combination of eq. ([3]) is fully affected by the node precessions 
of degree higher than 2, i.e. £ = 4,6,8,...; instead, the LAGEOS-LAGEOS II-LARES 
combination of eq. (|5]) will be fully impacted by the even zonals of degree higher than 4, i.e. 
£ = 6, 8, .... A realistic evaluation of the systematic uncertainty induced by the mismodeling 
in such uncanceled even zonals on the predicted Lense-Thirring signals, i.e. 



5f 



(aL)! 



< 



E 



gfii^) 



dJf 



5Ji,q = 2,3, 



(7) 



See also iRies et all fl2003f ): iPavlisI ( 120021 ) 
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Table 3: Nominal values of the coefficients ci of the present LAGEOS-LAGEOS II combina- 
tion, and ki, /c2 of the future LAGEOS-LAGEOS II-LARES combination according to Table 
[U The combined Lense-Thirring node precessions are 47.8 mas yr~^ (LAGEOS-LAGEOS 
II, eg. (El)), and 50.8 mas yr"^ (LAGEOS-LAGEOS II-LARES, eq. (IS))), respectively. 



Cl 



h 



k. 



0.5439211320 0.3603291106 0.0751007658 



has been the subject of several recent st udies summ arized in lloriol (120091 ). Concerning the 
present-day LAGEOS-LAGE OS II tests (lloriol 120091) . the total accuracy may be up to 2 — 3 
times larger than claimed by Ciufolini et al.l ( 120091 ): in the case of the futur e tests involving 
LARES, both gravitational ( lloriol l2009l ) and non-gravitational ( lloriol l2010l ) mismodeled 
perturbations should likely impact the mission at a level larger than the claimed 1%. 

In this paper we want to deal with another, subtle issue pertaining the systematic bias 
induced by the even zonal harmonics of the geopotentiajj. Indeed, all the studies performed 
so far relied upon the assumption of a perfect cancelation of J2 by the combinations of eq. 
([3]) and eq. (|5]). Actually, it depends on the accuracy with which their coefficients ci,ki, k2 
can be known; given the huge magnitude of the nominal J2-induced precessions of Table [2] 
with respect to the gravitomagnetic ones of Table [1], it has to be quite high to really allow 
for a measurement with a given level of uncertainty X%. Instead, until now, the coefficients 
of the combinations of eq. ([3]) and eq. (|5]) have always been computed with a few decimal 
digits. In other words, one has to evaluate 



5f 



(2L)| 



<5ci 



OLAGEOS II 
"J2 



and 



5f 



(3L)| 



I fcijfcj 



<5ki 



Q 



LAGEOS II 
J2 



+ 6k2 



qLARES 
J2 



(9) 



as further sources of systematic uncertainty with respect to the combined Lense-Thirring 
precessions which have to be added to Sf^'^^^j , g = 2,3. In order to avoid possible 
confusions and misunderstandings, it should be clarified that it would be incorrect to 
evaluate the impact of the uncertainties in the combinations' coefficients by only taking 
terms proportional to cross products of the errors like Sci6Ji, 6kiSJi, Sk2SJ£ instead of 
those proportional to 6ci,6ki,6k2 themselves, as done by us. Indeed, it is well known from 
elementary theory of errors that if an empirically determined quantity / depends on several 
parameters pj,j = 1,2... affected by uncertainties 6pj, the total uncertainty in / is just 



^/<E 



df 



dpj 



5pj 



(10) 



^I thank an anonymous referee of a previous paper of mine for having pointed out this 
issue to me. 
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Mixed terms of the form 

appear only in case of a correlation, which is absent in the present case. Indeed, the 
coefficients Ci,ki,k2 of the combinations f^^\q = 2,3 and the even zonals Ji of the 
geopotential are not solved-for parameterql^, simultaneously estimated in the same global 
solution: otherwise, one may look at their mutual correlations in the covariance matrix. 
Anyway, even if it was the case, a conservative evaluation of the total uncertainty would 
require to neglect the covariance by only retaining the linear sum of the individual 
mismodelled terms. 

The paper is organized as follows. In Section |2] we will deal with the ongoing 
LAGEOS-LAGEOS II tests. The LAGEOS-LAGEOS II-LARES case will be tackled in 
Section [3l while Section H] contains the summary and the conclusions. 



2. The LAGEOS-LAGEOS II case 

The coefficient ci of eq. (jlj) actually depends on the semimajor axes, the eccentricities 
and the inclinations of both LAGEOS and LAGEOS II. Thus, the accuracy with which 
it is possible to know it is set by the uncertainties in such Keplerian orbital elements. 
They are not directly measurable quantities being, instead, computed from the satellite's 
state vectors r and v whose components are, in turn, estimated in a least-square sense by 
processing the diffe rences between the observed and calculated ranges at different times 



(ITaplev et al.ll2004i ). 



Let us, now, consider in detail how to assess the uncertainty in the semimajor axis a 
due to a key geodetic parameter, i.e. the Earth's gravitational parameter GM which must 
be assumed as known to pass from the state vector to the Keplerian orbital elements. For a 
Keplerian orbit the semimajor axis is given by 

where r and v are the satellite's geocentric distance and speed, respectively. Thus, the 
relative uncertainty in a due to GM is 



6a 
a 



^,, ~ \gm) gm"" ~ \gm) i^^ • ^ ^^ 



^''Actually, the Lense-Thirring effect itself has never been explicitly modelled and solved- 
for in all the analyses performed so far. 



Averaging over one orbital period P^ = 27r/n by means of 



it turns out that 



dt _ (1 - e^)3/2 
^~ 27r(l + ecos/)2^ 

'6GM\ 



(^«Iga./) 



GM J 



Since fJMc Carthv fc Petitll2004f ) 



6GM 
GM 



2.00702 X 10"^ 



(14) 



(15) 



(16) 



the average uncertainties in the semimajor axes of LAGEOS, LAGEOS II and LARES are 
of the order of 

(^olageos) < 2.5 cm, 



(^flLAGEOS ii) < 2.4 cm, 



(17) 



(^flLAREs) < 1.6 cm. 

An issue is that the Earth's gravitational parameter GM is estimated by processi ng lon g 
SLR data sets in which just LAGEOS and LAGEOS II play a fundamental role ( jDunn 
20031 ). Moreover, the gravitomagnetic field of the Earth has never been accounted for in the 
solutions yielding GM produced so far, so that a twofold source of a-priori "imprinting" 
of the Lense-Thirring itself is present in the values of the Earth's GM adopted. It would 
be necessary to use figures obtained without including data from SLR targets, especially 
LAGEOS and LAGEOS II, although they may be less accurate. 

Actually, the total, realistic uncertainty in a should be even larger because of r and 
V entering eq. f lT^ . Indeed, concerning the uncertainty in r, it includes thd^^l cm— level 
accuracy in the station— satellite range p and the uncertainty in the geocenter— station 
position -Rsta, of the order of about 1 — 2 cm (JLejba fc SchillakI 120091 ). Anyway, in future 
calculation we will use the values of eq. ( TTTl) . 

Concerning the LAGEOS-LAGEOS II combination of eq. ([3]) used for the present-day 
tests, the uncertainty in its coefficient ci can be conservatively evaluated as 



Sci < 



dci 


6ai^ + 


dai^ II 


5aL II + 




SIl + 





5h 



L II- 



^^Of course, it has to be intended in the root— mean— square sense; it is not the mere 
single— shot mm— accuracy. 
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If, together with eq. fITTj) for the uncertainties in the semimajor axes, we assume a 
reasonable and reahstic value for the uncertainties in the inclinations of LAGEOS and 
LAGEOS II, i.e@ 51 = 0.5 mas, eq. (dH]) yields 



5ci = 1.30 X 10" 
corresponding to a residual J2 bias (see Table [2]) 



6ci 



OLAGEOS II 

"J2 



10.8 mas yr" 



(19) 



(20) 



and, thus to a_percent uncertainty in the Lense-Thirring combined signal of 23 %. Instead, 
if we considecj (5/l = 30 yuas, 51^ u = 10 /ias claimed by fICiufolini et al.ll2009[ ) we have 



6ci = 7.9 X 10"^ 
yielding an uncanceled J2 signal (see Table [2]) 



6ci 



qLAGEOSII 

J2 



6.5 mas yr 



-1 



(21) 



(22) 



which corresponds to a percent uncertainty of 14%. 



These results show that the issue of the imperfect cancelation of the largest node 
precessions due to J2 cannot be neglected in the evaluation of the total error budget, 
especially because the previous figures have to be added to those accounting for the 
mismodeling in the other even zonal harmonics of higher degree fully impacting the 
combination of eq. ([3]) . 



3. The LAGEOS-LAGEOS II-LARES case 

The case of the LAGE OS-LAGEOS II-LARES combination of eq. ([S]) can be treated in 
a similar way. Iloriol (120 lOl ) preliminarily dealt with it by considering the impact of clr and 
Jlr only on k2. Instead, one has to fully take into account the uncertainties of the orbital 
elements of LAGEOS and LAGEOS II as well in both ki and k2 according to 



6ki/2 < Yl 



dk 



1/2 



dipi 



Sipi, ipi = ol, /l, cll II, II II, olr. Ilk- 



(23) 



^^Indeed, it corresponds to a reasonable 6r ^ 1 + 2 = 3 cm from 61 ~ 6r/a. 

^^Such figures seem to be unrealistic because they would imply an accuracy 6r ^ a6I in 
reconstructing the orbits of LAGEOS and LAGEOS II, on average, of 0.2 cm and 0.06 cm, 
respectively. 
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By using eq. f iTTj) for 8a and assuming 51 = 0.5 mas for LAGEOS, LAGEOS II and LARES, 
the uncertainties in ki and k2 are 



6ki = 1.1 X 10' 



which yield an uncanceled J2 signal 



6ki 



(SLAGEOS II 
"J2 



+ 6k2 



Sko 



n LARES 
"J2 



2 X 10' 



(24) 



14.7 mas yr '^. 



(25) 



It corresponds to 29% of the combination of the Lense-Thirring node pre cessions. If, 



instead, in addition to eq. ( !T7|1 one adopts 5Jl = 30 /xas, (5/l ii = 10 /ias (ICiufolini et al. 
2OO9I ) and, say, (5/lr = 20 /ias the uncertainties in ki and /c2 are 



6ki = 5x 10"^ 6k2 = 2x 10"^. 
They yield a residual J2 signature 

Ski 



OLAGEOS II 
"J2 



Sk, 



qLARES 

"J2 



7.9 mas yr \ 



(26) 



(27) 



amounting to 16% of the predicted Lense-Thirring trend. 

Concerning the imperfectly canceled J4 signal, it turns out that it is of no concern 
amounting to 0.006 — 0.008 mas yr~^. 

Thus, independently of the lingering uncertainty in how to realistically assess the bias 
due to the mismodeling in the uncanceled even zonal harmonics of higher degree impacting 
in full the combination of eq. ([5]), the imperfect removal of the effect of J2 alone is sufficient 
to make dubious the achievement of the goal of a ~ 1% total accuracy in the future 
LAGEOS-LAGEOS II-LARES tests. 



4. Summary and conclusions 

One of the major sources of systematic uncertainty in the measurement of the 
gravitomagnetic Lense-Thirring precessions of the nodes Q of the laser-tracked LAGEOS- 
type satellites in the gravitational field of the Earth is given by the much larger competing 
classical node precessions induced by the even zonal harmonic coefficients Ji,i = 2,4,6, ... 
of the expansion in multipoles of the non-spherically symmetric terrestrial gravitational 
potential. The strategy followed so far to partially circumvent such an issue consisted of 
suitably designing linear combinations of the nodes of more than one satellite to cancel 
out, by construction, the effects of J2, as in the ongoing LAGEOS-LAGEOS II test, and 
of J4 as well, as in the future LAGEOS-LAGEOS II-LARES scenario. In addition to the 
usual systematic uncertainty due to the mismodeling in the even zonals of higher degree 
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which fully impact such combinations, another source of non-negligible uncertainty of 
gravitational origin has to be taken into account. It is due to the imperfect cancelation of 
the effects of J2 because of the uncertainty in the coefficients entering the combinations set 
up just to remove it. Indeed, the numerical values of such coefficients, released with just a 
few decimal digits so far, explicitly depend on the numerical values of the semimajor axes a, 
the inclinations / and the eccentricities e of the satellites involved. Thus, the uncertainties 
with which such Keplerian orbital elements are known unavoidably have repercussions onto 
the coefficients themselves. For uncertainties in the semimajor axes of 1 — 2 cm and of 
about 0.5 — 0.01/0.03 milliarcseconds in the inclinations we have shown that the resulting 
systematic bias due to the imperfect removal of the J2 signal may be as large as 14 — 29% 
of the Lense-Thirring signatures. 
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